Cylindrical Korteweg-de Vries solitons on a 
ferrofluid surface 



Dirk Rannacher, Andreas Engel 

Universitat Oldenburg, Institut fur Theoretische Physik, 26111 Oldenburg, Germany 
E-mail: rannacherOtheorie . physik . uni-oldenburg . de 

Abstract. Linear and non-linear surface waves on a ferrofluid cylinder surrounding 
a current-carrying wire are investigated. Suppressing the Rayleigh-Plateau instability 
of the fluid column by the magnetic field of a sufficiently large current in the wire 
axis-symmetric surface deformations are shown to propagate without dispersion in the 
long wavelength limit. Using multiple scale perturbation theory the weakly non- linear 
regime may be described by a Korteweg-de Vries equation with coefficients depending 
on the magnetic field strength. For different values for the current in the wire hence 
different solutions such as hump or hole solitons may be generated. The possibility to 
observe these structures in experiments is also elucidated. 



PACS numbers: 75.50.Mm, 47.20.Ma 



Submitted to: New J. Phys. 



Cylindrical Korteweg-de Vries solitons on a ferrofluid surface 



2 



1. Introduction 

Solitons are among the most interesting structures in nature. Being configurations of 
continuous fields they retain their localized shape even after interactions and collisions. 
Observed originally long ago as stable moving humps in shallow water channels pQ 
they have been established since then in various physical systems including optical 
waveguides, crystal lattices, Josephson junctions, plasmas and spiral galaxies (for an 
introduction see Long lasting efforts to theoretically describe their intriguing 

properties have culminated in the development of the inverse scattering technique [3] 
which is among the most powerful methods to obtain exact solutions of nonlinear partial 
differential equations 

Particularly popular examples for solitons in hydrodynamic systems are the 
solutions of the Korteweg-de Vries equation (KdV) 

d t u(x, t) + 6u(x, t)d x u{x, t) + dlu(x, t) = , (1) 

with x standing for a space coordinate and t denoting time. With u representing the 
surface elevation of a liquid in a shallow duct this equation can be derived perturbatively 
from the Euler equation for the motion of an incompressible and inviscid fluid Ej. 
The one-soliton solution of $1} is given by 

u(x,t) = ^se C h 2 (^(x-ct)^ (2) 

which for all values of c > describes a hump of invariable shape moving to the right 
with velocity c. The amplitude of the hump is given by c/2 whereas L = Ijyfc is a 
measure of its width. 

A decisive prerequisite to derive (JIJ is that to linear order in the field u the system 
under consideration admits travelling waves u ~ e *( fc:c -<^) w yfo dispersion relation 

uj = Co k + 0(k 3 ) for k^O, (3) 

where Cq denotes the phase velocity. Intuitively the invariant shape of the soliton solution 
may then be understood as the consequence of a delicate balance between nonlinearity 
and dispersion at higher orders of both u and k j2] • 

In the present paper we investigate cylindrical solitons of KdV-type on the surface 
of a ferrofluid in the magnetic field of a current-conducting wire. In order to conserve the 
radial symmetry of the problem we neglect gravity. A possible experimental realization 
of to this situation is to surround the ferrofluid column with a non-magnetic fluid of the 
same density. In this case the hydrodynamics of this fluid has to be treated as well. 

Ferrofluids are stable suspensions of ferromagnetic nano-particles in Newtonian 
liquids and behave superparamagnetically in external magnetic fields |15j . In the 
standard setup of a horizontal layer of ferrofluid subject to a homogeneous magnetic 
field an additional term proportional to k 2 shows up in the dispersion relation (JHJ) JUj 
which inhibits the derivation of a KdV equation in this geometry. On the other hand, for 
a ferrofluid cylinder in the magnetic field of a current-carrying wire the magnetic force 
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may replace gravity and allows for dispersion free surface waves in the long wavelength 
limit. This in turn paves the way to derive a KdV equation for axis-symmetric surface 
deformations on the ferrofluid cylinder [TT| I12j. 

Due to surface tension a long fluid cylinder is unstable to surface modulations 
resulting eventually in disconnected drops (Rayleigh-Plateau instability). Before 
embarking on the study of travelling waves on the fluid surface therefore means have to 
be found to suppress this instability. Fortunately, this can also be accomplished with 
the help of the magnetic field [THj . 

An accurate experimental investigation of solitons in hydrodynamic systems is 
notoriously difficult due to the ubiquitous presence of dissipation. Most quantitative 
studies have been devoted to hump solitons in shallow channels of water [TEIIH] whereas 
recently also the detection of hole solitons on the surface of mercury have been reported 
[T£] . In our present setup either hole or hump solitary waves are possible depending on 
the value of the applied current. We therefore hope that the present theoretical work 
may also stimulate new experimental investigations. 

The paper is organized as follows. In section 2 we collect the basic equations and 
boundary conditions. Section 3 is devoted to a linear stability analysis of a cylinder 
of ferrofluid in the magnetic field of a current-carrying wire. Here we demonstrate 
the possibility to suppress the Rayleigh-Plateau instability and establish the dispersion 
relation (j^J) for axis-symmetric surface waves. In section 4 we derive the KdV equation 
by multiple scale perturbation theory with details of the calculation relegated to two 
appendices. Section 5 provides the explicit form of the one- and two-soliton solution and 
gives some estimates for possible experimental realizations. Finally, section 6 contains 
some conclusions. 



2. Basic equations 



We consider a cylindrical column of ferrofluid surrounding a straight, thin, long, current- 
carrying wire under zero gravity. The ferrofluid is modelled as an incompressible, inviscid 
liquid of density p and constant magnetic susceptibility \ surrounded by a vacuum. 
Although we will eventually be interested in the non-linear evolution of the surface 
profile of the fluid the assumption of a linear magnetization law M = %H is quite 
reasonable for experimentally relevant parameters as will be discussed in section 5. We 
use cylindrical coordinates (r, 9, z), with the z-axis pointing along the wire (see figure^). 
The magnetic field is given by 

H = 7T~ e e > (4) 

2nr 

where J denotes the current through the wire. Due to the field a magnetization M 
builds up in the ferrofluid. The corresponding magnetic force, F m = /j (MV)H attracts 
the ferrofluid radially inward. The equilibrium free surface of the ferrofluid is hence 
cylindrical with the radius denoted by R. Deviations from this shape are parametrized 
by a function ((z, 9, t) according to r = R + C{z, 9, t). 
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Figure 1. Schematic plot of the system under consideration. A current-carrying 
wire is surrounded by a ferrofiuid column with magnetic susceptibility x an< i density 
p (region 0) under zero gravity. Region (2) is a nonmagnetic medium of negligible 
density treated as vacuum. The dynamics of the deflection £ of the surface from the 
perfect cylindrical shape with radius R is the central quantity of interest. The vector 
n denotes the normal on the free interface R + £(z, 9, t). 



The velocity field v(r, 9, z) inside the ferrofiuid is determined by the continuity 
equation 

V-v = 0, (5) 

and by the Euler equation 

p<9 t v + p(vV)v = -VP + /i (MV)H. (6) 

Here P(r, 6, z) denotes the pressure. We will only consider situations in which the flow 
of the fluid is irrotational, 

V x v = . (7) 
It is convenient then to introduce a scalar potential for the velocity 

v = V$ (8) 
which due to J5J) fulfills the Laplace equation 

A$ = (9) 
The Euler equation may now be integrated once to yield the Bernoulli equation 

(10) 



P — 2 ■ D ^ H 2 = const. 



p«9 t $ + |(V$) 2 + P 9 

The magnetic field has to obey the magnetostatic Maxwell equations [To] 
V-H = 
V x H = 0, 

both inside and outside the ferrofiuid. Denoting the respective fields by Hi and H 2 
equations (fTTj) allow the representations 

Hi = -Wi and H 2 = -W 2 (12) 

with the scalar magnetic potentials and ^2 also fulfilling the Laplace equation: 

A\&! = and A^ 2 = . (13) 



Cylindrical Korteweg-de Vries solitons on a ferrofluid surface 5 

Equations (jHJ), (jlUp . and (|13|) are to be complemented by boundary conditions. On 
the hydrodynamic side we have, assuming no radial extension of the wire, 

limd r $ = 0. (14) 

i — >o 

Moreover at the free surface we need to fulfill the kinematic condition 

d t C + d z $d z ( + = «9 r $ (15) 

as well as the pressure equilibrium [15] 

P = P + aK-^M 2 . (16) 

Here a is the surface tension, K := V ■ n denotes the curvature of the free surface, 
and M n is the magnetization perpendicular to the surface. The normal vector n on the 
surface is given by 

|V(r-C(M,*))l' 1 ' 

Note that ( = yields K = 1/R as it should be for the undisturbed cylinder. 
The boundary conditions for the magnetic field assume the form ^3] 



limd r ^ 1 = 



;is) 



lim d r ^ 2 = . 

r— *oo 

At the free surface we have 

n- V(* 2 - + = 

Equations (fT§|) describe the feedback of the flow of the ferrofluid onto the magnetic field. 

It is convenient to introduce dimensionless units. We measure all lengths in units 
of the cylinder radius R, and use the replacements 



(19) 



t^J&t, P^P, »^». (20) 

V a W p R 2,71 

The overall magnetic field strength which can be externally controlled by changing the 
current J is then characterized by the dimensionless magnetic Bond number 

Bo:=^. (21) 

Using the Bernoulli equation (fTUj) the pressure equilibrium (fT6j) at the free surface 
r = 1 + C(z, 9, t) is given by 

d t $ + ~ (V$) 2 + V • n - ^ ( X (n ■ V^) 2 + (V^) 2 ) = 1 - ^ . (22) 

Here the reference pressure Pq in (|16j) has been chosen such that $ = 0,^ = is a 
solution of 
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3. Linear stability analysis 

In this section we study the linear stability of the cylindrical interface given by 
£ = 0, $ = , = ^2 = 0. To this end we introduce small perturbations 
((9, z,t),4>(r,9, z,t),i[)i(r,0, z,t), and ^(r, 9, z, t) of the surface profile, velocity 
potential and magnetic potentials respectively and linearize the basic equations and 
their boundary conditions in these perturbations. From the translational invariance 
along the z-axis and equations (jOJ) and (fT3|) together with the boundary conditions (fT4*|) 
and (fTHj) it follows that these perturbations are of the form 

((9,z,t) = C n exp(in9 + ikz +pt) (23) 

<j>(r,9,z,t) = D n I n (kr) exp(in9 + ikz + pt) (24) 

ipi(r, 9, z, t) = A n I n (kr) exp(in9 + ikz + pt) (25) 

ip2( r , 9, z, t) = B n K n (kr) exp(in9 + ikz + pt) . (26) 

Here k denotes the wave number in z-direction, n 6 Z characterizes the azimuthal 
modulations, and p is the growth rate. The A n ,B n , C n and D n are constants (with 
their dependence on k and p suppressed) and I n (k) and K n (k) denote modified Bessel 
functions of order n |TJ|. 

Using the linearization of we may express A n and B n in terms of C n according 

to 

A = ■ KM r 

n mx i n {k)K' n {k) - ^m)Kn(k) n 

R= . W) {27) 

H mX I n (k)K' n (k) - Mk)Kn(k) ° n ' 

where the prime denotes differentiation with respect to the argument. In addition the 
linearized version of (j!5)) gives 

c -=wfk c - (28) 

Finally, linearizing (|22j) we find 

d t <P - dlQ - d 2 z C + (Bo - 1)C + Bodofa = (29) 
which when combined with (127)1 and ()28|) yields the dispersion relation 

rfW = {!-»'- Bo- e) + Ui J k f° — . (30) 

The reference state of a cylindrical column becomes unstable if combinations of k, n and 
Bo exist for which p n is positive. 

For Bo = we find back the well-known Rayleigh-Plateau instability accomplished 
by radially symmetric modes with n = 0. Modes with higher values of n are not able 
to destabilize the fluid cylinder. 

Since one has for all k 

m >0 and uk)K(k) 

lM I' n (k)K n (k) < ° (31) 
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we infer from (J30|) that p^(k) is a monotonically decreasing function of the magnetic 
Bond number Bo. The magnetic field hence always stabilizes the cylindrical surface. 
Consequently it may change the qualitative behaviour of the system only due to its 
influence on the n = modes. For n = the dispersion relation reads 

h(k) 



pl(k)=k^{l-Bo-k") . 
It is displayed for several values of Bo in figure (J2J). From 



(32) 

we see that the Rayleigh- 
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Figure 2. Plot of the square of the growth rate of axis symmetric distortions 
as a function of the wave number k as given by (|32[) for magnetic Bond numbers 
Bo = 0, 0.3, 0.7, 1 and Bo — 1.3 (from top to bottom). The inset shows a magnification 
of the region around k = 0. For Bo > 1 one has p 2 < for all k and no instability 
occurs. 

Plateau instability for a ferrofluid column will be suppressed by a sufficiently strong 
magnetic field fulfilling Bo > 1 [13J. Using typical parameter values as \ = 1.2, 
a = 0.03 J/m 2 and a fluid radius of R = 1 cm the system remains stable if the current 
exceeds the threshold J c ~ 89 A. 

It is instructive to investigate the dispersion relation (J30|) in the long wavelength 
limit (k 1). Using the expansion of the modified Bessel functions for small arguments 



[TH] we get 



Io(k) 



Hk)K{k) 

I' (k)K (k) 
Uk)K(k) 
I' n {k)K n {k) 



16 

^ if n > 



k 2 log k 
-1 if 



n > 



and hence find 



Pl(k) 

Pl(k) 



1 — Bo 



9- Bo 



16 



n(l 



n 



Bo) 



n 



k 4 + 0(k 

2 xBo 
X + 2 



0{k 2 ) if n>0. 



(33) 
(34) 
(35) 
(36) 

(37) 
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Therefore for n = and Bo > 1 the system exhibits surface waves ({z,t) ~ 
expi(kz — ut) with dispersion relation 

w(fc) = ^I*( 1 _^*.) + 0(*.) (39, 

The important point for what follows is that these surface waves become dispersion free, 
uj = Cofc, in the long wavelength limit k — > 0. The phase velocity is given by 

I Bo - 1 , , 

c o = V^— • (40) 
The situation is hence analogous to the shallow water equations which form the starting 
point for the derivation of the Korteweg-de Vries equation in a rectangular duct [Oj. 
Note that no such waves are possible for n > 1, cf. (|3*H|). 



4. Korteweg-de Vries equation 

In the previous section we have seen that the system under consideration admits to 
linear order in the surface deflection ( cylindrical, axis symmetric surface waves with 
no dispersion in the long wavelength limit k — > 0. From the experience with plane 
surface waves on shallow water jHj it is hence tempting to investigate whether at higher 
orders in k and ( nonlinear waves may be obtained for which the effects of nonlinearity 
and dispersion exactly balance each other. This could then give rise to axis symmetric 
soliton solutions in the present cylindrical geometry. 

In this section we show that it is indeed possible to derive a KdV equation for the 
surface deflection ((z, t) ^1] by using a multiple scale perturbation theory similar to the 
case of rectangular geometry. To this end we first observe that for an axis symmetric 
free surface the magnetic field problem decouples from the hydrodynamics and we have 
the exact result 

#i = 0. (41) 

This in turn implies V#i = (0, 1/r, 0) and therefore n ■ V#i = 0. Using moreover the 
explicit expression for n in terms of ((z,t) resulting from (j!7|) we get from 



ft* + \ [(^) 2 + + J - \j^- 2 = 1 - f . (42) 

The kinematic condition (|T5|) simplifies to 

d t ( + d z $d z ( = «9 r $ (43) 

The KdV equation appears in the limit of small nonlinearity, ( < 1, and small 
dispersion, fc< 1 with the proper balance between these two ingredients occurring for 
( = 0(k 2 ). To make this combined limit explicit we introduce a small parameter, e, and 
use the rescalings 

z t 
z^^, r^r, C^<, * -> 7=, $ -> Vtc $ , (44) 
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where Co is defined by (|4U|) . To derive the KdV equation we will need the two basic 
equations (jj2)) and up to order e 2 . Plugging (J4*4"j) into these equations we find to 
the required order 

d t $ + ~ [e (<9,$) 2 + («9 r $) 2 ] + 2C - -^(3Bo - 2) C 2 - - 2 d 2 z C = (45) 

Z ZCg Cg 

and 

e^C + e 2 9 2 $^C = «9 r $. (46) 

To get a suitable expansion for the velocity potential $ we note that from the 
Laplace equation © and the boundary condition (fPfj) one may derive the following 
representation for $(r, z,t) (see Appendix A) 

$(r, M) = £ r 2m e m ±-L^ %™%(z, t) (47) 

m=0 ^ *' 

with the so far undetermined function $> (z,t). 

Using this expansion for $ in ()45|) and (j4*o^) and observing that both equations hold 
at the interface, i.e. for r = 1 + e(, we get to the desired order in e 

e aa 2 A e r a ^^2 , e(35o-2) 



and 

1 



d t $ + 2C = - 2 (^ $ o) 2 + V 2g2 ; C 2 + T2 d% (48) 



<9*C + ~^$o = -ed t $ d z ( - -C d 2 z % + -«9 Z 4 $ (49) 
z z lb 



It is convenient to differentiate (|48|) with respect to z and to introduce the z-component 
of the velocity of the ferrofluid u = d z §. We then find the final set of equations to 
determine ( and u 

d t u + 2d z ( = e (\d t d 2 z u - ud z u + 3jB °~ 2 (d z ( + J (50) 

\4 Cn Cn / 



and 



<9*C + = e f- U cU - d z u + ^c\ 3 M ) . (51) 



(52) 



We now solve these equations perturbatively using the ansatze 

C{z, t, t) = ( (z, t, t) + eCi(z, t, r) + 0(e 2 ) 

u(z, t, t) = u (z, t, t) + eu x (z, t, r) + 0(e 2 ) , 

where we have introduced a second, slow time variable r := et. Plugging these 
expansions into ()50)1 and (|5l"| we find to zeroth order in e 

m \ / 



Co / V o 

where the linear operator L is given by 

d t 2d z 



2 



d z d t 



(53) 



(54) 
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The solution are dispersion free travelling waves of d'Alembert form 
u (z,t,r) = 2f(z-t,r) 

(55) 

(o{Z,t,T) = f(z-t,T), 

with a so far unspecified function f(x,r) where we have restricted ourselves to waves 
travelling to the right. 

To order e we find 
L fu x \ _ / -d T u + \d t d 2 z u Q - u d z u + ^pCocUo + jzd 3 z (o \ 
\ Ci / \ -drCo ~ u d z ( - ±( d z u + ±d 3 z u J 

This inhomogeneous equation involves again the linear operator L which is singular, 
cf. (|53|). Hence the inhomogeneity of this equation must be orthogonal to the zero 
eigenspace of the adjoint operator L + . The determination of L + and the projection of 
the r.h.s. of ()56j) onto the eigenfunction of L + with eigenvalue zero is done in appendix 
B. The solvability condition for (jHfij) finally acquires the form 

W + ^/ W + ^/ = (57) 

Using ([55)1 . denoting £o simply by £ and reversing the scalings (|44|) then yields the 
following KdV equation for the surface deflection ((z,t) 

dtC + c d z ( + 2 -^p^ (d z ( + dlC = (58) 

4co ozco 

When discussing the implications of this equations one has to keep in mind that it is 
valid for small ( only. 

5. Results 

Equation (|58|l is of the form 

dtC + c d z ( + Cl (d z ( + c 2 dlC = , (59) 
with the coefficients 



Bo -I 25o-3 , Bo -9 

Co = V^' Cl = ^^' and C2 = ^oV m 

all depending on the magnetic field strength Bo. From section III we know that we must 
have Bo > 1 since otherwise the fluid cylinder is susceptible to the Rayleigh-Plateau 
instability. Hence both c\ and c 2 may change sign for allowed values of Bo. 
The one-soliton solution of ()59)1 is of the form (cf. ((%))) 



C (D M ) = ^ sech 2 U^(z- (c + c )t)j . (61) 

where c < 1 is a free constant having the same sign as Ci. For Bo < 9 we have 
hence c < and the soliton has a slightly smaller velocity than the linear waves. If 
Bo < 3/2 also c\ < and therefore the amplitude of the soliton is positive, i.e. we have 
a hump soliton as shown in figure Et- For 3/2 < Bo < 9 on the other hand C\ > 
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and consequently (|UTj) describes a depression or hole soliton as depicted in figure Eb- 
Finally, for Bo > 9 we have C2 > 0, hence c > 0, and also Cj > 0. The soliton 
amplitude is therefore positive again and its velocity is now slightly larger than that of 
the corresponding linear waves. 



Figure 3. Schematic plot of a cylindrical hump-soliton (a) and hole-soliton (b). The 
parameter values are specified in the main text. 

To get some impression of the accessibility of the solution in experiments the 
results for the following parameter sets may be helpful. For a ferrofluid with x — 1-2, 
p = 1.12 g/cm 3 and a = 0.03 J/m 2 forming a cylinder of radius R = 1cm a current 
I = 100 A corresponds to Bo ~ 1.27. A soliton with amplitude A = 2 mm has then a 
velocity ofU = 1.8 cm/s and the width of the hump is about L = 20 cm. This soliton will 
hence be difficult to observe in an experiment. For a current of I = 294 A corresponding 
to Bo ~ 11 the extension reduces for the same amplitude to L — 1.6 cm with the 
velocity increasing to U = 12.3 cm/s. A hole-soliton with amplitude A = —2 mm, 
velocity U = 8.4 cm/s, and width L — 2.1 cm can be realized with a current of I = 235 A 
corresponding to Bo ~ 7. The latter two solitons are shown schematically in figure El 
Both should be easily observable experimentally. Note that for such values of the current 
the magnitude of error assuming \ = const, of the magnetization is about 5 %. 

A two-soliton solution may be derived using, e.g., Hirotha's method [Tl]. Depending 
on the value of the magnetic bond number one may combine either two hump or two 
hole solitons. The case of two hump solitons is described by the solution 

(V(z,t) = 8 V7 7 2 (62) 



where 



for i = 1, 2. 



7i = ^— an d 6 = exp ( A / — (z - z 0i - (q + c )t) ) (63) 




A snapshot of the solution is displayed in figure 01 its time evolution is characterized 
by figure El The main feature is the passing of the slower soliton by the faster one. After 
the interaction process the two solitons reemerge undisturbed which is the defining 
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property of a soliton solution. An animated version of the two-soliton solution is shown 
in a movie for parameters values as in figure 0] 




Figure 4. Schematic snapshot of the two-soliton solution described by eq. I|tj2fl for 
the parameter values R = 1.8 cm, Bo = 11, A\ = 8 mm, = 4 mm at t = —0.5. 




Figure 5. Evolution of the two-soliton solution ea. (|62J) at with time t = 
{ — 1.5, —0.5,0.5, 1.5,2.5}. The parameters are the same as in figure^] 

6. Conclusion 

In the present paper we have investigated nonlinear waves on the cylindrical surface of a 
ferrofluid surrounding a current- carrying wire under zero gravity. We have shown that 
for a sufficiently large current a Korteweg-de Vries equation for axis-symmetric surface 
distortion can be derived. Accordingly the system under consideration is well suited to 
experimentally investigate cylindrical solitons of KdV type. 

In order to observe these solitons first of all the ubiquitous Rayleigh-Plateau 
instability has to be suppressed. This can be accomplished by the magnetic field if 
the current exceeds a critical value which for experimentally relevant parameters is of 
about 100 A. 

We have shown that non axis-symmetric perturbations of the surface always 
disperse whereas axis-symmetric ones propagate almost dispersion-free if the wavelength 
is very large. Using the fact that for axis-symmetric surface deflections the magnetic 
field problem decouples completely from the hydrodynamic part a Korteweg-de Vries 
equation can be derived. The parameters in this equation depend on the magnetic 
field strength which gives rise to qualitatively different soliton solutions like hump and 
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hole solutions for different values of the current in the wire. The one- and two-soliton 
solutions were discussed in detail and conditions for their experimental realization were 
given. 

It should be noted that several approximations were used in our theoretical 
analysis. First of all the derivation of the KdV-equation is perturbative and therefore 
approximate as is typical for the theoretical discussion of solitons in hydrodynamic 
systems. Furthermore we have neglected the hydrodynamic influence of a non-magnetic 
fluid surrounding the magnetic column which is necessary to ensure zero gravity 
in experiments. Finally, viscosity was neglected throughout the analysis since the 
Korteweg-de Vries equation results from the Euler equation describing inviscid fluids. 
In experiments one will hence always see a damping of the soliton solutions with time 
due to dissipation by viscous shear flow [El[T3[THj. It is an attractive idea to counter- 
balance these viscous losses by appropriately chosen time dependent magnetic fields, 
however we were not able to find a suitable geometry for this idea to become operative. 
In any case a theoretical analysis aiming at this goal has to go beyond the quasi-static 
version of ferro-hydrodynamics employed in the present analysis and has to include 
magneto-dissipative couplings, see e.g. (2D]- 

We finally note that our system is an experimentally accessible realization of the 
introductory example for a soliton given in chapter 1.4 of the book by Lamb [21]. There 
an incompressible fluid inside a cylinder made of independent elastic rings is considered. 
The rings are supposed to deform axis-symmetrically in reaction to the fluid pressure. 
However, although confining the liquid tightly they must be uncoupled in order not to 
sustain elastic waves by themselves. Gravity is neglected altogether. Using the conser- 
vation of mass and momentum of the fluid and linear elasticity for the rings it is then 
possible to derive a KdV equation for axis-symmetric deformations of the rings. As we 
have shown in the present paper the somewhat unrealistic properties of the elastic rings 
can be mimicked by a cylindrical magnetic field if the fluid to be confined is a ferrofluid. 
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Appendix A. Expansion of $(r, z,t) 

After the rescalings the Laplace equation for the velocity potential $ takes the form 





r 



Representing $ as a power series in r 




(A.2) 



in 
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r m [(m + 2) 2 $ m+2 + e «9 2 $ m ] = 

m 

leading to the recursion relation 

$771+2 = -e- 



(m + 2) 2 

Because of the boundary condition (JT3jl we have 
^mr m " 1 $ m = 



14 



(A.3) 



(A.4) 



(A.5) 



implying $i = 0. From (|A.4|) we hence find $2m+i = for all m. The velocity potential 
may therefore be expressed entirely in terms of $o an d its derivatives 



00 f_i "i™ 



m=0 



d> (z, t) - 6 ^d 2 z $ Q (z, t) + ^di%{z, t) + 0(e 3 ) . (A.7) 



which coincides with f|4*7jl . 



(A.6) 



Appendix B. The solvability condition 

Under the usual scalar product 



(\&|*) = lim 



Z i>T 

dz I dt ^* ■ * 



(B.l) 



z,t^oo 4zt y_ z y_ T 

with \I/ = (it, C) we find for L + 

-ft -|5 
-2ft -ft 

The complete eigenmode to zero eigenvalue of L + is hence given by 

u (z,t) = f(z-t) -g(z + t) 

( (z,t) = 2(f(z-t)+g(z + t)), 

where / and g are arbitrary functions of a single argument. Setting the projection of 
the r.h.s. of on this mode equal to zero we find 

1 



(B.2) 



(B-3) 







lim 

Z,T^oo AZT 



Z pT 

dz I dt 



2ft/ + ^fdj - ^-Adlf )[f-g 



+ 2( -d T f-3fd z f + -d 3 z f)[f + g 



lim 



1 



, dz dt 

Z,T^oc AZT J_ z J_ T 



+ 



ABo-3 



/ft/ 



Wo- 11 



ft 3 / b 



• (B.4) 



Cylindrical Korteweg-de Vries solitons on a ferrofluid surface 
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The part of the integrand involving g may be written in the form g(z + t)d z F(z — t,r). 
Substituting £ = z — t,r] — z + t one realizes that these terms do not contribute for 
Z, T — > oo. Since moreover / is an arbitrary function of its argument (|B.4J) implies ()57jl . 
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